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Abstract 

Classical Molecular Dynamics (MD) simulations for a one-component plasma 
(OCP) are presented. Quantum effects are included in the form of the Kelbg 
potential. Results for the dynamical structure factor are compared with the 
Vlasov and RPA (random phase approximation) theories. The influence of 
the coupling parameter T, degeneracy parameter pA 3 and the form of the pair 
interaction on the optical plasmon dispersion is investigated. An improved 
analytical approximation for the dispersion of Langmuir waves is presented. 
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I. INTRODUCTION 



The model of a classical one-component plasma has - due to its simplicity - been widely 
investigated both theoretically and with various numerical and simulation methods, see e.g. 
} t]-[3|] and MM, respectively. Since the pioneering numerical work of Brush, Sahlin, Teller 

the thermodynamic and dynamic characteristics of the classical OCP have been stud- 
ied in detail. In particular, the dependence of the properties on the coupling parameter 
T = Aire 2 1 '(f/c^T), where f = (jf-) 1 ^ 3 is the mean interparticle distance and p the density, 
have been investigated up to very large values of T pW] . Among the most important ther- 
modynamic results is the prediction of crystallization at values of V of the order of 172-180 
PJlO[| . Furthermore, investigations of the dynamic properties of strongly correlated classical 
plasmas have indicated that the wave number dependent plasmon dispersion changes from 
monotonic growth, common for weakly coupled plasmas, to a decreasing dispersion around 

On the other hand, there is growing interest in the dynamic properties of dense quantum 
plasmas, in particular in astrophysics, laser plasmas and condensed matter. While the case 
of strong degeneracy (strong quantum limit) and weak coupling at very high densities is 
well described by the random phase approximation (RPA, see e.g. [|TT| -|T4}| ) , the properties 
at intermediate coupling and degeneracy remain poorly explored. Especially, one is interested 
in the dynamic plasma behavior in cases where the average kinetic energy is of the same order 
as the mean potential energy, i.e. T ~ 1, where collisionless theories such as the RPA fail, 
e.g. |i2l , |Il^ , |T5| . For these situations, quantum molecular dynamics (QMD) simulations [16 



are the appropriate numerical approach which, however, is yet lacking the required efficiency 
For weakly degenerate plasmas, with pA 3 < 1, where A is the DeBroglie wave length (see 
below), it is expected that one can perform much simpler classical MD simulations using 
effective quantum pair potentials, e.g. |p],pT7|| . These potentials can be derived from the 2- 
particle Slater sum using Morita's method. It is the aim of this paper to explore this MD 
approach in detail, especially for the analysis of the optical (Langmuir) plasmon dispersion. 
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It is natural to start this analysis with OCP-simulations because they have the advan- 
tage of the absence of a collapse of oppositely charged particles at small distances. On the 
other hand, the existence of a homogeneous background of oppositely charged particles leads 
to some additional technical difficulties compared to 2-component systems, due to restricted 
carrier rearrangement causing less effective screening of the Coulomb interaction. One major 
problem of MD simulations of dynamical properties is that the behavior at small wave num- 
bers is difficult to investigate. The reason is that large box-sizes are required which, for the 
analysis of high density plasmas, translates into large particle numbers. The current increase 
of available computer power gives one the possibility to investigate size-dependent properties 
like the density-density correlations (pg(O)pg(t)) for smaller k- vectors than before. In this 
paper, we are able to present accurate results for the dynamical properties of the OCP, such 
as the dynamical structure factor and the wave vector dispersion of Langmuir oscillations. 
Our simulations for intermediate values of the coupling parameter, T = 1 ... 4, show an in- 
teresting dispersion: the frequency increases up to a maximum and, for large wave numbers, 
decreases again. Further, we investigate the role of quantum effects by comparing simula- 
tions with the Coulomb potential and an effective quantum pair potential (Kelbg potential 
FT] ]) for the region of small and intermediate coupling. We found that quantum diffraction 



effects have noticable influence on the behavior of the optical dispersion curves. Increase of 
the degeneracy leads to a softenig of the dispersion o>(/c), especially at intermediate wave 
vectors. 

II. DYNAMICAL PROPERTIES OF THE OCP 

A. Statistical approach 

A central quantity to determine the dynamic properties of charged many-particle systems 
is the frequency-dependent dielectric function e(k,u) which, for the OCP, is given by 

e(k,u) = l-Uc(k)U(k,u). (1) 
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Here Uc(k) is the spatial Fourier transform of the Coulomb potential, Uc(k) = 4ire 2 /k 2 , and 

— # 

TL(k, uj) is the longitudinal polarization function. Thus, many-body effects enter the dielectric 
function via EL There exist many approximations for the latter function, the simplest one 
being mean-field theories which neglect short-range correlation effects, i.e. collisions between 
the particles. For the classical OCP, the mean-field result is the Vlasov polarization: 

„ ,r* . If d 3 v k dF(v) 

Here 5 — > +0, indicating the retarded (causal) character of the polarization and the dielectric 
function. Further, F is the distribution function. The Vlasov polarization applies only to 
classical plasmas, where the wave character of the particles can be neglected. Quantum 
effects are important if the interparticle distance or the Debye radius become comparable to 



the DeBroglie wave length A = h/ ^/2wmk~BT . Therefore, quantum diffraction effects should 
show up in the dielectric properties at large wave numbers. The quantum generalization of 
the Vlasov polarization is the RPA polarization function given by 

In this paper we consider only plasmas in equilibrium, so F and / are the Maxwell and 
Fermi function, respectively. One readily confirms that, in the limit of long wavelengths, 
k — > 0, indeed the RPA result (|3]) goes over to the Vlasov polarization function (Q). An 
important quantity which follows from dielectric function (P via the fluctuation-dissipation 
theorem is the dynamical structure factor S(k, uj) 

S(k,u;) = -^—lm-^— (4) 

7iU c {k)uj e (k,u) 

which shows the frequency spectrum of density fluctuations for a given value of k. 

As mentioned above, the mean field expressions (fj) and (|3]) neglect short-range correla- 
tions and are, therefore, valid only for weakly coupled plasmas, r < 1. There exist many 
theoretical concepts to go beyond the RPA which are based on quantum kinetic theory, 
density functional theory and other approaches. This is beyond the scope of this paper, see 



e.g. Ref. [ I2|, 131J15] and references therein. Here, we consider the alternative approach to 
the OCP at finite coupling which is based on molecular dynamics simulations. 



B. Molecular dynamics approach to the dynamical properties 

The dielectric and dynamical properties of an interacting many-particle system are easily 
accessible from the density-density correlation function which is defined as 

A&,t) = ^(Pi(t)P-s(0)), (5) 
where iV is the number of particles, ptit) is the Fourier component of the density, 

N 

= (6) 
i=l 

which is computed from the trajectories r*j(t) of all particles. The dynamical structure factor 
is just the Fourier transform of the density- density correlation function 



S(k,u) = — I dte iuJt A(k,t). (7) 

Z7T J 

— oo 

Equation (|7|) can be directly compared to formula (|J) and, thus, allows for a comparison of 
the simulation results with the statistical theories. Furthermore, Eq. (|7p allows to investigate 
the influence of quantum effects on the dynamical properties and plasmon dispersion of 
the OCP. Variations of the interaction potential (see below) directly affect the particle 
trajectories and, via Eqs. (H)-(0), the dynamical structure factor. 

III. DETAILS OF THE MD-SIMULATIONS 

The simulations have been performed in a cube of length L containing N interacting 
electrons on a uniform positive background. For this system, we solved Newton's equations 
of motion containing all pair interactions which are derived from a total potential U(r), 
see below. As an algorithm of motion we used a second-order scheme in form of the Swope 



algorithm Since our simulations are performed in the microcanocial ensemble, the mean 



kinetic energy may change. Therefore, to maintain the chosen value of temperature and T, 
we applied scaling (renormalization) of all velocities at every second step. 

A central goal of our simulations was to study the influence of quantum effects. We, 
therefore, performed several simulations which used either a Coulomb potential or an effec- 
tive quantum pair potential (see below). To permit flexibility in the use of the potential, 
U was divided into a short-range and a long-range part, U = U sr + U h , where quantum 
effects influence only U sr , whereas the behavior at large distances, U lT , is dominated by the 
long-range Coulomb interaction. Let us first describe the treatment of the long-range term. 



A. Long-range interaction 

The long-range interaction was computed in standard way using periodic boundary con- 
ditions and the Ewald summation procedure p0| , pl| . As a result, the long-range potential 



is given by the Coulomb interaction in the main box and all image cells: 

N 

U h (r) =4vre 2 £VW ald (r%), (8) 

„ ^ n ^- l eii C [^\{r + nL)/L\] n - n - n ^ 5 ' n ^ 27 exp(- 7m 2 ) cos(2™r7L) 1 , n , 
^EwaldW- ^ ^— ^ + X, ^ (9) 

n=0 1 1 n=f=0 

where erfc is the complementary error function, L - the side length of the simulation cell 
and n - a vector of integer numbers which labels the periodic images of the simulation 
box. In this expression, the first term corresponds to a potential of particles with Gaussian 
broadened charge distribution around the electrons with a width of y/n, the second one - 
the compensating Gaussian distributions, and the last one accounts for the influence of the 
homogeneous background. It turns out that the second term in (Q) can be reduced to 2 loops 
(one over the particles and one over the vectors n in the reciprocal space) and is not very time 
consuming. The more complicated part is the first term which contains three loops. In case 
of a two-component plasma, a proper choice of the width of the Gaussian distribution and use 
of periodic boundary conditions greatly simplifies this term due to cancellations. In contrast, 
for an OCP, the background cancels the interactions only partially, "statically". As a result, 
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convergence of the sum is worse, and one needs to take into account all first neighboring 
image cells (total of 26) at every time step. The contribution of all neighboring cells except 
for the main one (0 < |n| < \/3) was computed, before the start of the simulations and 
stored in 3-dimensional tables for the potential and forces. During the simulations, we used 
3D-bilinear interpolation at every step to obtain the values of the potential and forces for 
intermediate distances. We found that 100 grid points in every direction are adequate, so 
the total size of the table was 10 6 elements. The particle interactions inside the main (n = 0) 
cell were evaluated directly at every time step without minimum image convention. 



B. Short-range interaction. Quantum effects 

Let us now discuss the short-range potential. As has been shown by Kelbg and co-workers 
|I7| , [T9"|| , quantum effects can be treated efficiently by an effective pair potential, the Kelbg 
potential: 

U KELBG (r,T) = 4^(l^E2tlim + ^ erfc(r/A) ) (10 ) 

where A = As a consequence of quantum effects, this potential differs from the Coulomb 
potential at small distances r < A and is finite at r = 0. Further, it is temperature-dependent 
via the thermal DeBroglie wavelength. The Kelbg potential can be regarded as the proper 
quantum pair potential following from the two-particle Slater sum S2 without exchange 
effects: 

lnS 2 = -^f^ + 0(T*). (11) 

It treats quantum diffraction effects exactly, up first order in T. Frequently other quantum 
pair potentials have been used, including the Deutsch potential |[22|| , which has the same 



value at r = but differs from the Kelbg potential at intermediate distances. As was 
mentioned by Hansen P3 |, symmetry effects do not have a big influence on the dynamical 



properties (although they give a major contribution to the static properties, especially for 
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the light mass components). Using the Kelbg potential (|10D, we can immediately separate 
the short-range part of the interaction, 

V(r, T) = W ( " 6Xp( ; r2/A2) + ^erfc(r/A)) , (12) 

which has been calculated together with the first sum of Eq. (g) using the interpolation 
table. 

The Kelbg potential contains just the lowest order quantum corrections (lowest order in 
e 2 ) and is, thus, accurate at small coupling, T < 1. Nevertheless, we expect that it correctly 
reproduces the influence of quantum effects also at intermediate coupling, T < 5. Further 
improvements are straightforward, e.g. by including exchange effects or by evaluating the 
full two-particle Slater sum. We note that the described numerical procedure applies to such 
improved quantum pair potentials as well, even if they are not given analytically. 



C. Thermodynamic and dynamical quantities 

Solving Newton's equations with forces derived from the total potential U sr + U h , we 
computed thermodynamic and static quantities, such as total energy and pair distribution 
function in usual manner. The results will be presented in the next section. Here we discuss 
some details on computation of the dynamical properties, as they require much more effort 
and computation time in order to achieve sufficient accuracy. 

To obtain useful results for the dynamical structure factor, requires simulation results 
in a sufficiently broad range of wave numbers and frequencies. Natural units of the wave 
number and frequency are 1/f and the plasma frequency u p i = \J Ane 2 p/m, respectively, 
which will be used in the following. The minimum wave number k m i n depends on the size 
L of the simulation box and thus, for a given density or coupling parameter, on the number 
of particles N. One readily verifies that k min = 2n/L = 27r(p/iV) 1 / 3 or, using dimensionless 
wave numbers, q min = k min f = (67r 2 / TV) 1 / 3 . Clearly, to reduce k min requires an essential 
increase of the number of particles in the simulation. 
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The simulation accuracy can be further increased by taking advantage of the isotropy 
of the plasma in wave vector space. Indeed, in equilibrium, the density-density correlation 
function and dynamical structure factor should only depend on the absolute value of the 
wave vector. On the other hand, the simulations yield slightly different results for different 
directions of the wave vector. Averaging over all results corresponding to the same absolute 
value of k allows to reduce the statistical error. For example, the minimum wave number 
kmin corresponds to directions of k along either the x-, y- or z-axis, cf. Eq. @, so we can use 
the average of the three. The next larger value is \f2k min) corresponding to the diagonals 
in the x-y, x-z and y-z planes. The third value, \/3k m i n , corresponds to the space diagonal 
and is not degenerate; consequently it carries the largest statistical error. This is the main 
reason for the fluctuations of the numerical results for the wave vector dispersion, see for 
example Fig. 5. 

Finally, to resolve the collective plasma oscillations, the duration of the simulations has to 
be much larger than the plasma period. Also, increased simulation times leads to a reduction 
of the noise. We found that times of the order of 250 plasma periods are adequate. 

IV. NUMERICAL RESULTS 

We have performed a series of simulations for varying values of T and pA 3 , using the 
Coulomb and Kelbg potential. Also, time step and particle number have been varied until a 
satisfactory compromise between accuracy and simulation efficiency has been achieved. The 
parameters of the runs chosen for the figures below are summarized in Table I. We mention 
that kinetic energy conservation in all runs (if velocity scaling was turned off) did not exceed 
0.1%. Also, the results for the total energy (not shown), in case of the Coulomb potential, 
agree very well with data from the literature. 

We first consider the pair distribution function g(r) for varying interaction potentials 
and parameter values. Fig. 1 shows g(r) for three values of the coupling parameter, T = 
0.5, 1,4. As expected, the Coulomb pair distribution function is close to the Debye-Hiickel 
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limit for small coupling, with increasing T, the deviations, especially around r = f, grow 
systematically. The Kelbg pair distributions practically coincide with the Coulomb functions 
for r > 0.6 f but deviate from the latter at small distances of the order of the thermal 
DeBroigle wave length A where quantum effects are important. Clearly, with increasing 
degeneracy, the ratio A/f increases, and the deviations extend to larger distances and grow 
in magnitude. With increasing T, the deviations become smaller since Coulomb effects 
dominate the behavior at small distances. 

Let us now turn to the dynamical properties. In case of an OCP, charge and mass fluc- 
tuations are identical because of the rigid opposite charge background. In our simulations, 
we have calculated the density- density correlation function (||) and, by numerical Fourier 
transformation, obtained the dynamical structure factor S(q, uj) for several (from 6 to 10, de- 
pending on the simulation) wave numbers, the values of which are determined by the size of 
the simulation box L (see above). The value of the smallest wave number is given in Table I. 
The frequency dependence of S(q, uj) for several wave vectors is presented in Figs. 2-4 for the 
Coulomb and Kelbg potentials. Also, the results of the mean-field models are shown. The 
peak of the structure factor is related to the optical plasmon (Langmuir mode) of the elec- 
trons, its position shows the plasmon frequency Q(k), its width - the damping of the mode. 
In the limit k 0, Q(k) u p i for all models. For increasing wave numbers, the width of 
the peak grows steadily, and it merges with the continuum of single-particle excitations, e.g. 
|13| , [14]] , therefore, no results for larger wave numbers are shown. 

Consider now the results for the plasmon dispersion more in detail, cf. Fig. 5. First, we 
discuss the mean field results (Q) which are calculated using the Vlasov and RPA polariza- 
tions, Eqs. (H) and (|3]), respectively. The Vlasov result was computed using the formulas 
given in the review of Kugler ||24|| , and for the RPA, a code was developed which accu- 
rately evaluates the pole integration in Eq. (|3]). Both approximations show the same general 
trend for small and intermediate wave numbers: with increasing wave number, the plasmon 
frequency and the damping increase. At large q, the dispersion exhibits a maximum and de- 
creases again. In all situations, the RPA yields a slightly smaller frequency than the Vlasov 
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result, whereas the damping values are very close to each other. 

Let us now turn to the simulation results. The Coulomb and Kelbg simulations have 
been performed for exactly the same parameters, except for N and run time (cf. Table I). 
(Notice that, in contrast to the Kelbg case, the Coulomb simulations depend only on T which 
can be achieved by various combinations of density and temperature). Comparison of the 
two simulations shows, cf. Fig. 2, that the results for the structure factors are very similar 
in case of small T. Peak positions and widths as well as the low and high frequency tails are 
very close to each other. The reason is obvious: since the potentials (and pair distributions, 
cf. Fig. 1) differ only at a small interparticle distances of the order of A, differences in the 
structure factor would show up only at k > 2n/A, which is about an order of magnitude 
larger than the wave numbers shown in Fig. 2. There, the plasmon peak has already a width 
of the order of the frequency and no longer describes a well-defined collective excitation. 

It is now interesting to compare the simulation results to the theoretical approximations. 
The first observation is that the simulation peaks are significantly broader, cf. Fig. 2. This 
is obvious since the simulations fully include interparticle correlations missing in the mean- 
field results. Consequently, the plasmon damping contains collisional damping in addition 
to the Landau damping (which is the only damping mechanism in the mean-field models). 
Correspondingly, the plasmon peaks in the simulations are shifted to smaller frequencies. 
This effect grows with increasing wave number as well as with increasing coupling (see also 
Fig. 5). We note that, in our simulations, this shift is observed for all wave numbers, which 



is in contrast to the result of Hansen [see Fig. 9 of Ref. [p3] for q = 0.6187]. In other words, 
the plasmon dispersion curves from the MD simulations are lower than the mean-field result 
for all wave vectors k, which is seen more clearly in Fig. 5. As expected the MD curves for 
the structure factor are much closer to the RPA than to the Vlasov result. 

In Fig. 5 we plot the optical plasmon dispersion curves for three values of the coupling 
parameter for the Vlasov and RPA dispersions together with the simulation results. We 
further show the well-known analytical approximation to the Langmuir dispersion, 
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/ 2X 1 / 2 

= Wp, ( 1 + V) • ( 13 ) 

Clearly, this predicts a monotonically increasing dispersion. However, this approximation is 
valid only for k < 1/rp and for V < 1. Let us now consider the simulation results which do 
not have this limitation. In Fig. 5 we show the MD results for a Coulomb potential and for the 
Kelbg potential for three values of the degeneracy parameter, pA 3 = 0.1, 0.5, 1.0. One clearly 
sees that, for these parameters, the dispersion is positive, duj(q)/dq > 0, up to wave numbers 
of the order of one over the mean interparticle distance. For larger q, the dispersion changes 
sign. This general trend is observed for the Coulomb and the Kelbg potential. On the other 
hand, with increasing quantum effects, pA 3 , the deviations between the two potentials are 
growing, which becomes more pronounced as T increases, cf. the curves for T = 1 and T = 4: 
the dispersion in case of the Kelbg potential shows a softer increase with increasing wave 
number and reaches a lower maximum approximately at the same wave number as in the 
Coulomb case. We mention that this sign change of the dispersion has not been reported 
by Hansen [[!]]. Comparing the simulations with the mean- field results, we again see that 
the MD dispersions proceed lower than the mean field results, and this effect grows with 
increasing T and increasing wave number. Once more, we confirm that the RPA dispersion is 
much closer to the MD result than the Vlasov dispersion, at least for T < 0.5. [As mentioned 
above, the simulation results for the dispersion show certain statistical fluctuations due to 
the varying accuracy of the results for the different wave numbers]. 

Let us now consider the plasmon damping more in detail. Fig. 6, shows the damping 
(full width at half maximum of the plasmon peak of the structure factor) as a function of 
wave number. It is interesting to compare with the familiar analytical expression from the 



Vlasov theory, e.g. [26 



ITT y/l + 3/t 2 13 

6 ( K ) = \h — ^ — e 2 ( 14 ) 



K 3 



where k = hrp is the dimensionless wave number in units of the inverse the Debye radius 
r£> given in Table I. Formula (|14]) is derived under the condition that the damping is much 
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smaller then the frequency [S(q) C and is limited to small wave numbers k < 1. As 

expected, the damping given by formula ( |i4|) which is only Landau damping, is much smaller 
than the damping found in the simulations, as the latter contain the full collisional damping 



also. Obviously, for small coupling and small q, Eq. ( |T4| ) shows the correct trend. However, 
deviations increase rapidly with growing coupling parameter. Furthermore, the simulations 
which are not limited to small wave numbers, show a qualitatively different behavior at 
large q: a monotonic increase of the damping. Interesingly, with increasing T the damping 
is reduced, cf. Figs. 6a,b. 

Finally, we try to extend the analytical result for the plasmon dispersion, Eq. (|13D, to 
larger T and to include quantum effects. To this end, we used the MD data with the Kelbg 
potential to construct an improved fit of the form oo(q) = u p i (1 + aq 2 + bq 4 ) 1 ^ 2 . The result 
is shown in Fig. 7 for T = 1 and T = 4. Due to the large fluctuations in the dispersion 
data and the increasing damping for large wave numbers, we used a weighted fit where the 
smallest q— values had the largest weight and the statistical errors of the individual points 
have been taken into account. Table II contains the resulting fit parameters. According to 
this data both parameters a and b depend on T and pA 3 . The parameter a is close to 1/r in 
agreement with Eq. (|T3D, but with increasing T, deviations are growing, compare Table II. 
We see no systematic influence of quantum effects on the parameter a for T = 1. Noticeable 
effects show up for T = 4, where increased degeneracy leads to a reduction of the coefficient 
a. The second fit parameter allows to qualitatively reproduce the change of the sign of the 
dispersion. The overall agreement is satisfactory for wave numbers up to the inverse mean 
interparticle distance up to which the plasmons are comparatively weakly damped. 

V. DISCUSSION 

We have presented classical molecular dynamics simulations of the dielectric properties 
of a one-component plasma at intermediate coupling and degeneracy, T < 4 and pA 3 < 1. 
While classical MD simulations can be extended to very large values of T, they have lim- 
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ited applicability to quantum plasmas. We used, as an effective quantum pair potential, 
the Kelbg potential which correctly describes quantum diffraction effects for small T. In 
general, we found that the simulation results for the dielectric properties and the plasmon 
dispersion with the Coulomb and the Kelbg potential are rather close, but start to deviate 
from each other as V increases. Nevertheless, the use of the Kelbg potential is preferable. It 
correctly treats close collisions, i.e. the two-particle interaction at distances smaller than the 
DeBroglie wavelength. This is of even higher importance in the case of two-component plas- 
mas where the Kelbg potential allows to avoid the collapse of oppositely charged particles. 
Therefore, the present investigation should be important for future work on two-component 
plasmas. Finally, we mention that the Kelbg potential is only the first term of a T expan- 
sion. Therefore, for T > 1 the account of higher order corrections to the quantum diffraction 
effects is necessary. Work on this subject is in progress. 
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TABLES 

TABLE I. Parameters of the molecular dynamics simulations with the Kelbg potential. Num- 
bers in parentheses refer the runs with Coulomb potential. 



L 


pA 3 


p, [cm -3 ] 


T, [K] 


u p i, [is] 1 


r D /r 


N 




run time, [T p i] 


0.5 


0.1 


9.12 TO 21 


1.126 TO 5 


5.387 


0.816 


500(250) 


0.491(0.619) 


515(341) 


0.5 


0.5 


2.28-10 23 


3.292-10 5 


26.940 




400(250) 


0.529(0.619) 


429(429) 


1.0 


0.1 


1.14-10 21 


2.228-10 4 


1.905 


0.577 


250 


0.619 


290(327) 


1.0 


0.5 


2.85-10 22 


8.23 TO 4 


9.524 




250 


0.619 


682(682) 


1.0 


1.0 


1.14-10 23 


1.31 TO 5 


19.048 




250 


0.619 


477 


4.0 


0.1 


1.78-10 19 


1.76-10 3 


0.238 


0.289 


250 


0.619 


570(227) 


4.0 


1.0 


1.78-10 21 


8.17 TO 3 


2.381 




250 


0.619 


716 



TABLE II. Fit parameters of the Langmuir dispersion curves shown on Fig. |^. The fit equation 
was taken in the form of uj{q)/u p i = (1 + aq 2 + bq 4 ) 1 / 2 . Parameters of the fit for L = 1 and pA 3 = 
0.1 are less reliable, because of the absence of data for big wave vectors, cf. Table I. 



L 


pA 3 


a 


b 


1.0 


0.1 


1.013 ± 0.031 


-0.260 ± 0.023 


1.0 


0.5 


1.074 ± 0.041 


-0.288 ± 0.013 


1.0 


1.0 


0.975 ± 0.055 


-0.259 ± 0.018 


4.0 


0.1 


0.169 ± 0.015 


-0.034 ± 0.006 


4.0 


1.0 


0.121 ± 0.007 


-0.025 ± 0.003 
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FIGURES 




FIG. 1. Pair distribution functions for T= 0.5 (upper figure), T= 1.0 (middle figure), 4.0 (lower 
figure), and pA 3 = 0.1, 0.5, 1.0 for systems with Coulomb and Kelbg potential. Further, the 
Debye-Hiickel (DH) limit is shown (solid line). Line styles are the same in all three figures. Inset 
in the middle Fig. shows the influence of the degeneracy at small distances. The result for T = 
4.0 , pA 3 = 0.1 with Kelbg potential are not distinguishable from the Coulomb result and are not 
plotted. 
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FIG. 2. Dynamical structure factor for an OCP at T = 1 and pA 3 = 0.1 from MD simulations 
with Coulomb and Kelbg potentials. In addition, Vlasov and RPA results are shown. The wave 
numbers are shown in the figures in units of f, i.e. q = kf. 
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FIG. 3. Same as Fig. 2, but for T = 0.5 and pA s = 0.1. The values of the wave numbers differ 
from Fig. 1 due to the different particle numbers, cf. Table I. 




FIG. 4. Same as Fig. 3, but for pA 3 = 0.5. 
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FIG. 5. Optical plasmon dispersion for various coupling and degeneracy parameters from MD 
simulations with Coulomb and quantum potentials. Also shown are results of the Vlasov and RPA 
approximations, and of the analytical approximation of Eq. (|l3|). For T = 4.0 and />A 3 = 0.1 
(upper graph) the MD simulations with Kelbg potential and the RPA curve are not shown since 
they almost coincide with the Coulomb simulation and the Vlasov curve, respectively. 
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FIG. 6. Damping of Langmuir waves from MD simulations with the Kelbg potential for various 
values of T and pA 3 . Solid lines are the analytical small damping limit of the Vlasov theory, 
Eq. (14). 
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FIG. 7. Dispersion of Langmuir oscillations from MD simulations with the Kelbg potential for 
various values of the coupling and degeneracy. Symbols are MD results, lines the best fits to the 
low wave number part (q < 1/f), the fit formula and parameters are given in Table II. 
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